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ABSTRACT 


Selective  Reenlistment  Bonuses  (SRBs)  are  offered  to 
improve  retention  in  designated  military  occupational 
specialties  (MOSs)  for  specified  years-of -service  intervals 
(zones).  The  amount  of  the  bonus  is  set  by  assigning  an  "SRB 
Multiplier"  for  each  MOS  and  zone  combination  (cell). 

Determination  of  multipliers  is  modeled  as  a  nonlinear 
knapsack  problem  which  is  then  linearized  to  a  generalized 
assignment  problem.  The  objective  is  to  minimize  the  sum 
over  all  cells  of  a  weighted  squared  deviation  from  the 
reenlistment  target  in  each  cell.  Lagrangian  relaxation 
provides  lower  bounds  and  feasible  solutions.  The  best 
feasible  solution  is  improved  using  a  greedy  heuristic  to 
apportion  unexpended  funds. 

A  FORTRAN  77  computer  program  implements  the  procedure. 
Data  for  FY86  yields  a  0- 1  integer  program  with  4895  binary 
variables  and  980  constraints.  A  solution  within  .01%  of 
optimality  is  obtained  on  an  IBM  3033AP  in  1.7  seconds  and 
on  an  IBM  PC  in  about  four  minutes. 
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I.  INTRODUCTION 


The  United  States  Marine  Corps  (USMC)  uses  Selective 
Reenlistment  Bonuses  (SRBs)  to  increase  retention  in  desig¬ 
nated  military  occupational  specialties  (MOSs)  and  years-of- 
service  intervals  (zones).  Those  MOSs  and  zones  which  will 
receive  bonuses  as  well  as  the  dollar  amounts  of  the  bonuses 
are  determined  by  the  assignment  of  "SRB  Multipliers."  In 
this  study,  the  problem  of  assigning  multipliers  given  a 
limited  budget  is  modeled  and  solved  as  a  mathematical 
program.  In  this  chapter,  the  SRB  program  is  described, 
terms  are  explained  and  the  problem  itself  is  defined.  In 
addition,  the  solution  approach  is  briefly  discussed  and  the 
structure  of  this  thesis  is  outlined. 

A.  PROBLEM  STATEMENT 

The  Marine  Corps  SRB  program  was  developed  as  a  reen¬ 
listment  incentive  to  increase  retention  in  designated  MOSs. 
In  effect,  an  SRB  is  a  sum  of  money  offered  to  individuals 
with  specific  skills  and  years  in  service  to  encourage  them 
to  reenlist  and  thereby  selectively  increase  manning  levels. 
The  bonus  is  applied  at  reenlistment  points  that  fall 
between  21  months  and  14  years  of  active  duty  service  (ADS). 
This  period  is  divided  into  three  reenlistment  zones.  As 
defined  in  the  applicable  Department  of  Defense  (DoD)  direc¬ 
tive  [Ref.  1],  the  period  between  21  months  and  6  years  of 
ADS  is  designated  as  Zone  A,  the  period  between  6  and  10 
years  as  Zone  B  and  the  period  between  10  and  14  years  as 
Zone  C. 

SRBs  are  not  offered  to  everyone  who  reenlists,  rather 
they  are  directed  at  selected  MOS  and  zone  combinations, 
called  "cells”  in  this  study.  For  each  cell,  there  exists  a 
reenlistment  target  which  is  defined  as  the  desired  number 
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of  reenlistments  in  that  cell.  In  many  cells,  enough  Marines 
will  reenlist  without  any  added  inducement  so  that  the  reen¬ 
listment  targets  will  be  automatically  met  or  exceeded.  In 
others,  however,  the  absence  of  a  bonus  will  result  in 
insufficient  reenlistments.  Through  the  offer  of  bonuses  to 
reenlist,  it  is  possible  to  reach  or  approach  the  targets. 
To  permit  additional  flexibility  in  awarding  bonuses,  the 
amount  of  an  SRB  can  vary  from  a  minimum  of  0  to  a  maximum 
of  $16,000.  From  experience,  it  is  known  that  larger 
bonuses  result  in  greater  retention  rates  although  the  rate 
at  which  retention  increases  with  the  size  of  the  bonus 
differs  greatly  between  cells.  The  total  amount  of  bonus 
money  that  can  be  awarded  per  fiscal  year  is  limited  by  the 
SRB  program's  annual  budget.  Since  an  individual  receives 
75%  of  the  bonus  as  a  lump  sum  payment  with  the  remainder 
apportioned  in  equal  annual  payments  over  the  term  of  the 
reenlistment,  the  total  SRB  budget  must  consist  of  two 
parts.  The  first  is  money  owed  from  bonuses  from  previous 
years  and  its  apportionment  is  non-discretionary .  The  rele¬ 
vant  budget  here  is  the  remainder  which  can  be  allocated  for 
current  year  bonuses . 

Once  the  decision  is  made  to  consider  the  use  of  a  non¬ 
zero  SRB  to  improve  manning  in  a  particular  cell,  the  actual 
monetary  size  of  the  bonus  is  based  on  increments  of  the 
eligible  member's  monthly  base  pay.  The  individual  receives 
an  amount  equal  to  his  monthly  base  pay  multiplied  by  that 
cell's  SRB  multiplier  and  the  number  of  years  of  additional 
obligated  service,  subject  to  a  maximum  of  $16,000.  The  DoD 
directive  stipulates  that  the  multipliers  may  not  exceed  6 
and  may  be  fractional  amounts.  The  Marine  Corps'  implementa¬ 
tion  of  this  policy  is  somewhat  more  restrictive.  It 
chooses  not  to  use  fractional  multipliers  and  does  not 
permit  a  multiplier  greater  than  5  for  reenlistments  in  Zone 
A,  4  in  Zone  B  and  3  in  Zone  C.  Assigning  a  multiplier  of  0 
is  equivalent  to  no  bonus. 


Presently,  SRB  multipliers  are  assigned  manually  in  an 
iterative  procedure.  Multipliers  are  selected  for  cells 
based  on  the  desired  number  of  reenlistments  for  those 
cells.  Once  a  particular  combination  of  SRB  multipliers  has 
been  selected,  a  microcomputer  based  spreadsheet  package  is 
used  to  determine  if  the  proposed  combination  is  feasible, 
i.e.,  within  the  allowed  budget.  If  it  is  not,  multipliers 
are  changed  and  the  calculation  is  repeated. 

At  this  time,  there  is  no  objective  function  with  which 
alternative  sets  of  proposed  multipliers  can  be  evaluated  or 
compared.  Instead,  the  officer  responsible  for  assigning 
multipliers  relies  exclusively  on  his  experience,  judgement 
and  various  rules  of  thumb  to  judge  the  des i rab i 1 i 1 t y  of  a 
particular  set  of  multipliers.  Even  if  a  sensible  objective 
function  were  defined,  it  is  doubtful  that  an  optimal  solu¬ 
tion  could  be  determined  using  a  manual  procedure.  Of  the 
approximately  1000  cells,  on  average,  450  are  assigned  non¬ 
zero  SRB  multipliers  during  any  given  period.  It  is  very 
likely  that  the  dimensionality  of  the  problem  would  quickly 
overwhelm  any  ability  to  manually  arrive  at  an  optimal 
allocation. 

This  study  proposes  an  alternative  procedure.  It  is 
automated  due  to  the  large  amount  of  data.  An  objective 
function  is  defined  and  justified,  thereby  permitting  the 
comparison  of  different  sets  of  multipliers.  Further,  an 
optimization  method  built  around  this  objective  function  and 
the  budget  constraint  is  described  and  tested.  Using  this,  a 
nearly  optimal  set  of  SRB  multipliers  can  be  determined. 
This  will  result  in  a  more  effective  allocation  of  the  SRB 
program  budget  . 

B.  PROBLEM  SCOPE 

The  goal  of  the  SRB  program  is  to  reduce  expected  short¬ 
falls  in  the  number  of  reenlistments  in  particular  cells  by 
increasing  the  retention  rate  in  those  cells  through  the 


offer  of  a  reenlistment  bonus.  The  model  that  is  developed 
in  this  study  tries  to  achieve  that  goal  by  determining  a 
set  of  SRB  multipliers  which,  by  a  selected  definition  ,  can 
be  considered  optimal.  This  is  a  single  period  model;  no 
attempt  is  made  to  forecast  sets  of  multipliers  for  subseq¬ 
uent  years  based  on  expected  retention  resulting  from  the 
current  year's  set  of  multipliers.  It  allocates  that  part 
of  the  budget  not  obligated  for  annual  payments  for  bonuses 
from  previous  years  and  does  not  consider  the  effect  of  the 
25%  of  the  bonus  that  is  apportioned  in  future  years. 
Furthermore,  the  model,  through  its  use  of  expected  values, 
is  essentially  deterministic.  Other  assumptions  regarding 
the  model  are  explicitly  addressed  in  the  succeeding 
chapter . 

C.  MODEL  AND  SOLUTION  APPROACH 

In  the  mathematical  formulation  that  is  proposed  to 
model  the  assignment  of  SRB  multipliers,  the  objective  func¬ 
tion  to  be  minimized  is  the  sum  of  individual  cell  penal¬ 
ties.  These  penalties  are  nonlinear  increasing  functions  in 
the  deviation  from  the  reenlistment  target  for  the  cell. 
Inclusion  of  the  budget  constraint,  which  is  also  nonlinear, 
results  in  a  nonlinear  knapsack  problem.  The  nonlinearity, 
coupled  with  other  complicating  factors,  prevents  the 
problem  from  being  efficiently  solved  using  common  knapsack 
algorithms.  For  descriptive  and  computational  purposes,  the 
problem  is  reformulated  as  a  special  generalized  assignment 
problem . 

In  this  generalized  assignment  formulation,  cells  repre¬ 
sent  the  tasks  that  must  be  assigned.  Unlike  conventional 
generalized  assignment  problems,  there  is  but  a  single  agent 
who  must  accomplish  all  the  tasks  but  who  has  alternative 
processes  for  each.  These  alternative  processes  correspond 
to  the  various  multiplier  levels  permitted  for  that  cell. 


Because  of  the  special  structure  of  this  problem, 
Lagrangian  relaxation  is  particularly  well  suited  as  a  solu¬ 
tion  technique.  Lagrangian  relaxation  has  been  successfully 
used  in  scheduling  (e.g.,  [Ref.  2]),  location 
(e.g.,  [Ref.  3])  and  set-covering  problems  (e.g.,  [Ref.  4]) 
as  well  as  generalized  assignment  problems  [Ref.  9],  Of 
interest  is  that  previous  applications  to  generalized 
assignment  problems  have  shared  the  characterist ic  of 
multiple  tasks  and  multiple  agents;  a  review  of  the  litera¬ 
ture  has  not  revealed  an  application  with  a  formulation 
similar  to  the  one  developed  in  this  study. 

An  important  advantage  to  the  Lagrangian  relaxation 
approach  is  that,  in  the  process  of  establishing  a  lower 
bound  on  the  optimal  solution,  feasible  solutions  are  also 
obtained.  A  heuristic  method  is  developed  which  improves 
the  best  feasible  solution  uncovered  in  the  bounding  proce¬ 
dure.  The  heuristically  improved  set  of  multipliers  is 
accepted  as  the  final  solution  provided  its  objective  value 
is  sufficiently  close  to  the  lower  bound. 

D.  THESI'  OUTLINE 

This  thesis  develops  and  presents  a  method  for  deter¬ 
mining  SRB  multipliers  in  the  USMC.  In  Chapter  II,  the 
problem  is  formulated  first  as  a  nonlinear  knapsack  problem. 
Then,  to  facilitate  the  description  of  the  solution  approach 
and  for  computational  purposes,  it  is  reformulated  as  a 
special  generalized  assignment  problem.  The  solution  meth¬ 
odology  is  presented  in  Chapter  III.  The  approach  uses  the 
technique  of  Lagrangian  relaxation  combined  with  a  heuristic 
procedure  to  provide  a  final  solution.  Details  on  the 
implementation  of  this  process  and  computational  results  are 
presented  in  Chapter  IV.  Conclusions  and  recommendations 
are  contained  in  Chapter  V.  Lastly,  a  listing  of  the  source 
code,  user  instructions  and  samples  from  input  and  output 
files  are  included  as  appendices. 


II .  MODEL  FORMULATION  AND  DESCRIPTION 

In  this  chapter,  the  process  of  determining  multipliers 
for  all  cells  is  formulated  first  as  a  simple  nonlinear 
integer  program.  Coefficients  are  defined  and  described  and 
the  rationale  behind  the  selection  of  the  objective  function 
is  explained.  For  explanatory  and  computational  purposes, 
the  problem  is  then  converted  to  a  special  case  of  the 
generalized  assignment  problem. 

A.  MODEL  DEVELOPMENT 

The  following  formulation  is  developed  to  model  the 
assignment  of  SRB  multipliers: 


Indices : 

i  =  1 ,  .  .  .  ,m 
Data : 


Functions 


ci(xi) 

ri ( xi ) 


Decision  Variables: 


Formulat ion : 


cells 

budget 

penalty  for  setting  multiplier  at  x^ 
in  cell  i 

cost  of  setting  multiplier  at  x^ 
in  cell  i 

maximum  permissible  multiplier  for 
cell  i.  n^  =  5,4,3  for  Zone  A,B,C 
respect ive ly 

mu ltiplier  for  cell  i 


i...; :  pried  to 
v  ;  .1*..  i uns  ip. 
st.  hers.  The 


rents  a  weighted  squared  devia- 
of  reenlistments  that  would  be 
: e r  x -  were  offered.  It  is  t h e 
The  first  term  repre- 

each  cell.  A  weighting  scheme 
some  colls  are  considered  more 


second  set  of  t e rm s 


Q  L  ( 


-D,  )*]-} 


:  s  *  he  i  red  .lev  i  it.  ion  from  the  desired  number  of  reen- 

;  ;  s  in  cell  but  where  Q  may  be  used  to  weight  the 

nt-na  I ;  y  i  r.curred  by  overages  ( the  second  difference  term) 
di  f  ft- rent  ly  f  r-mi  the  base  penalty  incurred  by  shortages  (the 
first,  difference  t  e  r:n )  .  The  Q,  the  "over/under"  factor,  has 
been  set  to  1  in  this  study  but  values  somewhat  less  than  1 
;r. ;  gbit  be  appropr  ivit.  e  . 

Squaring  the  deviation  has  the  effect  of  penalizing 
large  deviat ions  much  more  severely  than  smaller  ones.  This 
is  appealing  because  it  evens  out  shortages  and  overages 
among  cells.  Most  would  agree  that  it  is  less  detrimental 
to  the  force  as  a  whole  to  "spread  the  shortages  around" 
than  to  have  a  potentially  debilitating  shortfall  in  a  few 
cells.  Spreading  overages  around  is  desireable  for  logis- 
t.  i  ca  1  reasons  . 


The  deviation 

terms  are 

ca 1 cu  1  a t  ed 

by  squaring 

the 

d  1  f  f 

erence  between 

the  desire 

d  number  of 

reenl i stment s 

Di 

and 

t.  he  expected 

number  of 

reenlistments  R^(x^)E^. 

The 

f  unc 

t l on  Rx ( x  L  )  is 

disc  re  t  e 

and  defined 

only  at  integer 

values  of  on  the  interval  [0,n^].  Specifically,  when  the 
individuals  in  cell  i  are  offered  an  SRB  multiplier  x^,  a 
fraction  R^(x^)  of  the  Ej  individuals  eligible  to  reenlist 
will  choose  to  do  so.  The  values  of  R^(x^)  have  been  esti¬ 
mated  from  responses  to  previous  bonus  offers;  more  sophis¬ 
ticated  methods  for  determining  those  functional  values,  for 
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example,  by  including  the  effect  of  past  and  present 
economic  conditions,  are  currently  being  studied  [Ref.  5]. 

As  stated  previously,  otherwise  identical  deviations  in 
different  cells  are  not  viewed  equally.  Generally  it  costs 
the  Marine  Corps  considerably  more  to  replace  or  train  an 
individual  in  a  highly  technical  MOS  than  in  a  less  tech¬ 
nical  MOS .  Further,  the  actual  number  of  Marines  in  a  cell 
influences  the  acceptability  of  shortages.  For  example,  a 
shortfall  of  5  in  a  cell  containing  1000  could  be  more 
easily  tolerated  than  that  same  shortfall  in  a  cell 
containing  only  10. 

Therefore,  each  cell  i  is  assigned  the  weighting  coeffi¬ 
cient  is  the  "training  index"  for  cell  i  which 
is  just  the  cost  of  training  a  Marine  in  cell  i  normalized 
by  the  maxin.um  training  cost  over  all  cells.  A^~^  is  the 
reciprocal  of  the  actual  number  of  Marines  in  cell  i.  is 
an  exogenous  weighting  factor  which  permits  the  user  of  the 
model,  if  uesired,  to  exercise  his  judgement  in  weighting 
cells  or  blocks  of  cells;  presently,  all  Wi  are  set  to  1. 

The  form  of  the  resource  function  r^(x^)  is  straightfor¬ 
ward.  This  function  represents  the  expected  number  of 
dollars  that  must  be  allocated  in  order  to  offer  cell  i  an 
SRB  multiplier  of  x^  and  is  the  product  of  the  expected 
number  of  reenlistments  in  cell  i  that  would  result  from 
offering  a  bonus  x^,  R^(x^)E^,  and  the  size  of  the  bonus 
min  { ( . 75 )PiSixi ,  $12,000}  where  $12,000  is  751  (the 
percentage  of  the  bonus  that  the  individual  receives  as  a 
lump  sum  payment)  of  the  maximum  permissible  bonus  amount  of 
$16,000. 

Problem  PI  is  recognizable  as  a  nonlinear  knapsack 
problem  in  which  cells  represent  commodities  and  the  budget 
represents  the  weight  or  cube  constraint.  However,  the 
nonlinearity  in  the  objective  function  and  budget  constraint 
make  solution  by  standard  branch  and  bound  techniques 


impossible  [Ref.  6],  and  make  solution  by  dynamic  program¬ 
ming  inefficient  since  standard  reduction  techniques 
[Ref.  7]  are  inapplicable.  A  straightforward  dynamic 
programming  solution  is  possible  using  the  cells  as  stages 
and  dollars  as  the  sole  state  variable.  However,  signifi¬ 
cant  computational  difficulties  would  arise  due  to  the  large 
number  of  stages  (on  the  order  of  10  )  combined  with  a  state 
variable  that  for  FY86  can  achieve  any  integer  value  on  the 
interval  [0,7x10^].  The  solution  methodology  that  is  used 
could  be  applied  directly  to  the  nonlinear  knapsack  formula¬ 
tion.  However,  the  explanation  of  the  implementation  of  this 
methodology  will  be  more  transparent  if  the  problem  is  first 
converted  to  a  linear  integer  problem,  in  particular,  a 
generalized  assignment  problem.  The  following  section 
describes  this  new  formulation. 

B.  CONVERSION  TO  A  GENERALIZED  ASSIGNMENT  PROBLEM 

A  typical  formulation  for  the  generalized  assignment 


problem  is: 

[Ref. 

8] 

Indices : 

k  =  1,, 

.  .  .  ,n 

agent 

i  =  1,  • 

.  .  .  ,m 

task 

Data : 

cik 

cost  if  task  i  is  assigned  to  agent  k 

rik 

amount  of  resource  expended  if  task  i  is 

assigned  to  agent  k 

ak’bk 

minimum  and  maximum  amounts  of  the 

resource  that  may  be  expended  by  agent  k 

Decis ion 

Variable 

s  : 

1  if  task  i  is  assigned  to  agent  k 
0  otherwise 

Fo  rmu lat ion : 

min  XXcikxik  (p2> 

I  k 


a 

■ 


,  * 


*  4 


a 


a 


ft 


B 


a 


■ 


15 


l 


s  .  t .  a 


k  ^  Lrikxik  ^  bk 


Ixik  =  1 

k 


k  =  1 ,  .  .  .  ,n 


i  =  1 ,  in 


xik  e  {0,1} 

In  contrast  to  the  classical  assignment  problem,  the 
first  n  constraints  imply  that  more  than  one  task  i  can  be 
assigned  for  completion  to  agent  k  provided  that  the 
resource  constraints  on  agent  k  are  not  violated.  (Only  one 
type  of  resource  is  allowed.)  The  second  set  of  m 
constraints  requires  that  each  task  be  assigned  exactly 
once  . 

In  the  SRB  multiplier  selection  problem  the  model  is 
simplified  to  the  following: 


Indices : 

i  =  1, 
j  =  0, 

Data : 


cells 

SRB  multiplier  values 

ci( j  ) 
r±(j  ) 

budget 

maximum  permissible  multiplier  for  cell  i 
n-  =  5,4,3  for  Zone  A,B,C,  respectively 


Decision  Variables: 


1  if  multiplier  j  is  selected  for  cell  i 
0  otherwise 


Formulat ion: 


in  XLcijxij 


( P3a ) 


x-jj  e  {0,1} 


Here,  a  single  agent  is  required  to  perform  all  m  tasks  but 
has  at  his  disposal  n^  alternative  processes  for  each.  The 
above  formulation  is  completely  equivalent  to  PI  although  a 
continuous,  i.e.,  linear  programming  (LP),  relaxation  of 
this  formulation  might  be  very  weak  since  neither  r^(x^)  nor 
c^(x^)  will  typically  be  convex  functions.  No  problems 
arise,  however,  since  no  LP  relaxation  is  employed  in  the 
solution  procedure. 

This  formulation  of  the  model  is  quite  general.  The 
coefficients  c^j  and  r—  are  directly  calculated  by  evalu¬ 
ating  the  functions  c^(x^)  and  r^(x^)  °f  the  preceding 
section  at  all  feasible  values  of  x^.  Provided  that  the 
assumption  of  independence  between  cells  is  retained,  alter¬ 
native  functions  are  permissible.  It  will  be  seen  in  the 
succeeding  chapter  that  the  solution  methodology  is  equally 
general  and  introduction  of  alternative  functions  is 
trivial . 


III.  SOLUTION  METHODOLOGY 


With  the  problem  and  formulation  now  defined,  this 
chapter  details  the  approach  taken  to  find  a  satisfactory 
set  of  SRB  multipliers.  Commonly,  bounds  for  integer 
programs  are  established  using  an  LP  relaxation;  the  special 
structure  of  this  problem  and  the  likely  weakness  of  the  LP 
relaxation,  however,  suggests  an  alternative  approach  using 
Lagrangian  relaxation.  The  relaxed  formulation  is  extraor¬ 
dinarily  easy  to  solve  as  an  integer  program  and  the 
resulting  bound  is  at  least  as  good  as  that  produced  using 
an  LP  relaxation.  Furthermore,  in  the  course  of  determining 
the  bound,  feasible  sets  of  multipliers  are  encountered,  the 
best  of  which  is  heurist ically  improved  to  yield  a  final 
solution . 

A.  LAGRANGIAN  RELAXATION  TO  ESTABLISH  A  LOWER  BOUND 

The  technique  of  Lagrangian  relaxation  is  well  suited 
for  those  integer  programming  problems  that  would  be  rela¬ 
tively  easy  to  solve  were  it  not  for  one  or  more  compli¬ 
cating  constraints.  Moreover,  it  has  been  successfully  and 
efficiently  applied  previously  to  generalized  assignment 
problems  [Refs.  8,9]. 

As  described  in  Fisher's  excellent  primer  on  Lagrangian 
relaxation  [Ref.  10],  the  approach  is  to  move  complicating 
constraints  into  the  objective  function  using  the  product  of 
the  Lagrangian  multiplier  and  the  constraint  violation  as  a 
penalty  term.  To  demonstrate,  consider  the  linear  integer 
programming  problem: 

min  cx  (PM 
s  .  t  .  Ax  ^  b 


Dx  ^  e 

x^O  and  integer 


where  the  constraints  Ax^b  are  the  complicat ing  ones.  This 
is  relaxed  to: 


minx  cx  +  >v(Ax-b)  (P5a) 

s  .  t .  Dx  ^  e 

x^O  and  integer 

where  it  is  assumed  that  this  is  easy  to  solve  for  a  fixed 
row  vector  >»^0.  For  such  a  X,  an  optimal  solution  to 
problem  P5  provides  a  lower  bound  on  the  objective  value  for 
P4.  The  best  possible  lower  bound  from  a  relaxation  of  the 
form  of  P5  is  found  by  solving: 


rnax^^o  minx  cx  +  A.(Ax-b) 
<  s  .  t  .  Dx  ^  e 


x^O  and  integer 


(P5b) 


Furthermore,  this  bound  is  at  least  as  good  as  that  provided 
by  an  LP  relaxation  [Ref.  11]. 

Examination  of  the  generalized  assignment  formulation 
developed  in  the  preceding  chapter  indicates  that  the  only 
truly  complicating  constraint  is  the  budget  constraint; 
integer  optimization  with  this  constraint  relaxed  consists 
of  m  separable,  "multiple  choice"  problems  which  are 
trivial.  Thus,  the  relaxed  formulation  with  scalar  X  is: 


where  the  inner  portion  of  the  objective  function  may  be 
equivalently  written: 


min  VVx-  ,  (c,  ,+/.r-  J-A.B  (eqn  3.1) 


r  f  —r 


For  a  fixed  >»,  the  inner  minimisation  is  performed  by 
selecting  a  minimum  (c  — *>.r^j)  for  each  cell.  The  outer 
maximisation  problem  is  a  convex  optimisation  problem  which 
is  easily  solved,  since  >„  is  a  scalar  variable,  by  first 
bracketing  X  and  then  performing  a  bisection  search. 

The  bracket  about  the  optimal  X  is  readily  established. 
The  lower  bound  is  zero.  For  X  sufficiently  large,  the 
solution  to  the  inner  minimization  is  to  expend  as  few 
dollars  as  possible,  i.e.,  xi0=1  for  all  i.  L  ,  the  upper 
bound  on  interesting  values  of  ,  is  the  smallest  value  of  X 
for  which  this  solution  is  optimal.  Therefore,  LQ  is  the 
smallest  value  of  X  satisfying: 

ci0  +  ^ri0^cij+^-rij  V  (eqn  3.2) 

Recalling  that  r^j  represents  the  cost  of  offering  cell  i  an 
SRB  multiplier  j  and  that  j=0  corresponds  to  no  bonus,  r^g^O 
for  all  i.  Thus,  equation  3.2  becomes: 

ciOScij+^rij  V  ( eqn  3.3) 
and ,  as  a  result , 

Lu  =  max{  (ciQ -cij  )/rij  )  V  i,j^l  (eqn  3.4) 

For  >V  =  LU,  a  tie  exists  between  the  solution  x^g  =  l  for  all  i 
and  at  least  one  solution  where  x  — =  1  for  some  i  and  some 
j^l.  For  unusual  data,  the  latter  solution  could  be  infeas¬ 
ible.  To  ensure  that  at  least  one  feasible  solution  is 
obtained,  in  practice  Lu  is  replaced  by  (l+c)L  where  c  0. 

The  solution  to  P3b  provides  a  lower  bound  on  the 
optimal  solution  to  the  SRB  multiplier  problem.  Upper  bounds 
are  provided  by  feasible  sets  of  SRB  multipliers  encountered 
while  solving  P3b.  Such  solutions  arise  since,  for  X 


ft 


sufficiently  large,  the  entire  budget  will  not  be  consumed. 
The  best  of  the  feasible  sets  of  SRB  multipliers,  that  is, 
the  one  with  the  smallest  objective  function  value  and, 
correspondingly,  the  lowest  upper  bound  value,  is  improved 
using  the  heuristic  procedure  described  in  the  following 
section. 

B.  MARGINAL  RATE  OF  RETURN  HEURISTIC 

It  seems  reasonable  to  expect  that  the  best  feasible 
solution  obtained  in  the  process  of  optimizing  P3b  would 
require  the  expenditure  of  nearly  the  entire  SRB  budget. 
There  might,  however,  be  some  budget  dollars  remaining,  a 
residual  that  could  be  apportioned  to  further  reduce  the 
objective  value  of  the  best  feasible  solution  which  is  the 
present  upper  bound.  Accordingly,  a  heuristic  method  for 
allocating  this  residual  was  developed  using  the  concept  of 
marginal  rate  of  return. 

Given  a  best  feasible  set  of  SRB  multipliers  from  the 
Lagrangian  procedure,  the  incremental  cost  for  each  cell 
that  would  be  realized  by  increasing  that  cell's  multiplier 
by  1  is  calculated  for  all  those  cells  with  multipliers  not 
already  at  their  maximum  values.  If  this  incremental  amount 
is  less  than  the  unallocated  portion  of  the  budget  it  is 
possible  to  increase  that  cell's  multiplier  by  1  while 
remaining  feasible.  Let  such  a  cell  have  x^j=l.  Then,  the 
marginal  rate  of  return  for  the  cell  is  defined  to  be: 


RORi_  ^cij  _ci  ,j  +  l^  /  (ri,j*-l“rij) 


( eqn  3.5) 


This  quotient  expresses  the  improvement  realized  in  the 
objective  function  per  dollar  spent  when  cell  i's  SRB  multi¬ 
plier  is  increased.  In  the  heuristic  procedure,  the  cell 
with  the  maximum  positive  ROR^  has  its  multiplier  incre¬ 
mented.  The  residual  portion  of  the  budget  is  reduced  and 
the  process  is  repeated  until  it  is  no  longer  possible  to 


increase  Che  multiplier  of  any  cell.  The  resulting  set  of 
SRB  multipliers  is  accepted  as  the  final  solution  and  upper 
bound  for  the  SRB  multiplier  problem  provided  that  the  ratio 
between  the  lower  and  upper  bound  is  close  to  1. 

Other  heuristic  improvement  procedures  are  possible  but 
were  not  implemented  in  this  study.  For  example,  by  simul¬ 
taneously  incrementing  one  cell's  SRB  multiplier  by  1  and 

O 

decrementing  another's  (an  0(m‘1)  operation)  it  might  be 
possible  to  improve  on  the  upper  bound  while  remaining 


IV .  IMPLEMENTATION  AND  COMPUTATIONAL  RESULTS 

A  FORTRAN  77  computer  program  was  written  to  implement 
the  solution  methodology  described  in  the  preceding  chapter. 
Using  data  obtained  from  the  Manpower  Department  at 
Headquarters,  USMC  (HQMC),  the  model  was  run  on  both  an  IBM 
3033AP  and  an  IBM  PC.  To  judge  the  robustness  of  the  model, 
additional  testing  was  performed  using  randomly  perturbed 
data . 

A.  COMPUTER  PROGRAM 

The  program  that  performs  the  Lagrangian  bounding  proce¬ 
dure  and  which  yields  the  heur i s t i ca 1 ly  improved  final  solu¬ 
tion  was  written  in  ANSI  standard  FORTRAN  77.  No  functions 
unique  to  the  system  were  used.  When  the  program  was  run  on 
the  IBM  303 3AP  it  was  compiled  by  the  IBM  VS  FORTRAN 
compiler  at  0PT(2);  on  the  IBM  PC  it  was  compiled  using 
Ryan-McFar land ’ s  RM/FORTRAN  compiler  with  code  optimization. 
The  program  requires  six  input  files:  Parameter  file, 
Reenlistment  Plan  file.  Reenlistment  Rates  file,  Training 
Cost  file,  Cell  Size  file,  and  a  Cell  Weight  file.  File 
specifications  are  detailed  in  Appendix  D.  As  output,  the 
program  produces  a  file  containing  an  MOS/Zone  listing  with 
the  proposed  SRB  multipliers.  For  comparison  purposes,  the 
objective  value  as  well  as  the  percent  deviation  between  the 
final  solution's  objective  value  and  the  lower  bound  on  the 
optimal  objective  are  output. 

Under  some  circumstances  the  user  may  wish  to  preassign 
the  SRB  multipliers  in  specific  cells.  For  example,  the 
user  may  wish  to  declare  the  multipliers  for  all  three  zones 
in  MOS  2  112  to  be  some  specified  value.  This  may  be  accom¬ 
plished  by  "tagging"  those  cells  in  the  Cell  Weight  file  and 
entering  the  desired  SRB  multiplier  in  the  appropriate 


column  in  the  data  file.  Additional  details  are  provided  in 
Appendix  D.  Comparison  of  the  objective  values  before  and 
after  the  introduction  of  these  preassigned  multipliers 
enables  the  user  to  gauge  the  effect  of  a  prior:  multiplier 
ass ignment  s . 

B.  TEST  DATA 

The  FY86  data  that  was  used  to  develop  and  test  the 
model  was  obtained  from  the  Enlisted  Plans  Section  of  the 
Manpower  Department  at  HQMC .  Model  input  consists  of  those 
values  listed  in  equation  2.1  and  2.2  organized  into  five 
data  files  and  one  problem  parameter  file.  Because  the  rele¬ 
vant  data  is  extracted  from  a  variety  of  sources,  it  was 
deemed  desirable  to  place  them  in  different  files.  The 
information  contained  within  the  Parameter  file  is  primarily 
run  specific.  Samples  from  the  various  input  files  are 
contained  in  Appendix  B. 

C.  COMPUTATIONAL  RESULTS 

The  data  described  above  leads  to  a  problem  with  979 
cells.  This  in  turn  results  in  a  problem,  in  the  form  of 
P3b,  with  4895  binary  variables  and  980  constraints.  The 
objective  function  value  from  the  final  solution  is  within 
.01%  of  the  lower  bound  established  through  the  Lagrangian 
procedure.  In  the  process,  99.99%  of  the  budget  was 
expended.  Appendix  C  contains  a  sample  from  the  output 
file. 

Running  time  on  the  IBM  3033AP  was  1.7  seconds.  1.4  of 
the  1.7  seconds  is  used  to  read  in  the  data  and  calculate 
the  coef  ficients  c^j  and  r  —  .  The  Lagrangian  procedure  and 
the  heuristic  use  most  of  the  remaining  .3  seconds.  The 
time  required  to  write  the  solution  was  negligible. 

After  preliminary  development  and  testing  on  the  IBM 
3033AP,  the  program  and  data  were  copied  onto  a  floppy  disk 
ette  and,  using  the  RM/ FORTRAN  software,  installed  on  an  IBM 


PC  configured  with  512Kb  RAM,  an  8087  Math  Coprocessor  and 
two  360Kb  floppy  disk  drives.  5fo  changes  to  the  code  were 
required  in  moving  it  to  the  PC.  On  the  IBM  PC,  the  program 
runs  in  approximately  4  minutes.  Most  of  this  4  minutes  is 
consumed  reading  in  the  data  from  the  floppy  diskette; 
slightly  faster  times  might  be  expected  using  a  system  with 
a  hard  disk  arrangement. 

Additional  testing  was  conducted  to  verify  the  robust¬ 
ness  of  the  solution  procedure  and  to  confirm  that  the 
excellent  results  that  were  observed  were  not  merely  fortu¬ 
nate  happenstance  stemming  from  a  "good"  set  of  data.  In 
this  procedure,  the  data  read  by  the  model  was  perturbed  in 
this  manner: 


d ’ =  _d*UNIF ( . 5 , 1 . 5  )  j 
where 

d'  perturbed  datum 

d  original  datum 

UNIF(a,b)  random  variate  generated  from  the  continuous 

uniform  distribution  between  a  and  b 
x,  largest  integer  not  exceeding  x 

This  randomization  process  was  performed  on  all  the  data 
except  that  response  rates  were  capped  at  1  and  not  discre- 
t i zed . 

In  all,  50  additional  model  runs  were  conducted  in  this 
manner,  each  with  a  different  data  set.  Solution  times  on 
the  IBM  3033AP  varied  slightly,  ranging  between  1.5  and  1.9 
seconds.  In  each  case,  the  final  solution  was  within  . 02" 
of  the  lower  bound  on  the  optimal  solution.  Typically,  the 
heuristic  improves  the  best  feasible  solution  with  respect 
to  the  lower  bound  by  less  than  .04%.  In  a  situation  in  BR 

which  exceedingly  tight  bounds  are  not  necessary,  the 
heuristic  could,  in  fact,  be  eliminated. 


V  .  CONCLUSIONS  AND ^RECOMMENDATIONS 

In  this  thesis,  the  process  of  determining  Selective 
Reenlistment  Bonus  multipliers  subject  to  a  limited  budget 
was  modeled  as  a  mathematical  program.  In  the  course  of 
doing  this,  an  objective  function  based  on  expected  devia¬ 
tions  from  cell  reenlistment  targets  was  defined.  This  is  an 
important,  advance;  formerly  there  was  no  method  by  which 
alternative  sets  of  multipliers  could  be  compared.  Using 
this  objective  function,  a  procedure  was  developed  which 
determines  a  nearly  optimal  set  of  SRB  multipliers.  The 
speed  with  which  the  solution  is  found,  a  modest  core 
storage  requirement  and  a  compact  ,  specialised  FORTRAN  code 
proves  the  procedure  amenable  for  PC  use.  This  is  an  impor¬ 
tant  advantage  given  the  PC's  wide-spread  availability.  The 
inclusion  of  an  exogenous  weighting  factor  and  the  capa¬ 
bility  to  preset  multipliers  allows  the  user  additional 
flexibility  in  using  the  model. 

The  model  and  solution  method  that  are  employed  are  very 
general  and  only  require  that  the  assumption  of  inter- cell 
independence  be  maintained.  Other  penalty  and  resource 
utilization  functions  are  readily  incorporated  into  the 
model.  Improvements  in  the  estimates  for  the  response  rates 
should  further  increase  the  accuracy  with  which  the  model's 
multipliers  achieve  the  cell  targets. 

Provided  the  user  possesses  the  requisite  hardware  (the 
IBM  PC  configuration  explained  in  Chapter  IV  should  be 
considered  the  minimum)  and  a  FORTRAN  77  compiler,  the 
program  is  completely  operational.  A  copy  of  the  source 
code  is  listed  in  Appendix  A.  Other  pertinent  information, 
including  examples  of  input  and  output  files  and  user 
instructions,  is  contained  within  Appendices  B,  C,  and  D. 


The  future  portends  increasingly  restrictive  budgets 
within  the  Department  of  Defense  and  additional  pressure  to 
maximize  benefits  given  limited  resources.  With  respect  to 
maximizing  the  benefits  realized  from  the  SRB  program,  a 
nearly  optimal  assignment  of  SRB  multipliers  goes  far 
towards  achieving  that  goal. 


APPENDIX  A 

LISTING  OF  SOURCE  CODE 


IMPLICIT  REAL"  8  (A-H,0-Z) 

IMPLICIT  INTEGER"4  (I-N) 

REAL"8  C ( 1100 , 0 :  5)  ,A(  1100,0:  5) 

INTEGER*2  MOSfllOO) . IZONE ( I 100 ) , IXBEST ( 1 100 ) , IX ( 1 100 ' 
COMMON  CINF, BUDGET, TEST, NUMCEL 


CINF  = I . D22 

CALL  INITAL ( C , A , HLAM , MOS .IZONE . IX) 

CALL  BOUND (ZL , IX , ZU . IXBEST . HLAM . C . A , COST , MOS ) 
CALL  HEURISTC, A, IXBEST, ZH, COST, MOS) 

WRITE  (16,50) 


WRITE  (16,50) 

FORMAT  ( IX , ' ******************************** 

PCTDEV  =  100  .  DO '- ( 1 .  DO  -  ZL  /  ZH ) 

WRITE  (  16 . 100 )  ZH 
FORMAT  (IX. 'OBJECTIVE  FUNCTION  VALUE  IS:',F11.3) 
WRITE  (  16 ,  125J PCTDEV 

FORMAT  flX, 'THIS  SOLUTION  IS  WITHIN  ' ,F11.8, 

%  OF  THE  OPTIMAL') 

WRITE  (16,50) 

WRITE  (16,300) 

DO  150  1=1, NUMCEL 

WRITE  (16, 200 )MOS ( I ) , IZONE ( I ) , IXBEST ( I ) 
CONTINUE 

FORMAT  (IX, 15 , 3X , 12 , 9X , 12  ) 

FORMAT  (3X, 'm6s' ,3X, 'ZONE  ,5X, 'MULTIPLIER' ) 

STOP 
END 

SUBROUTINE  INITAL(C , A , HLAM , MOS , IZONE , IX ) 

IMPLICIT  REAL"  8  (A-H.O-Z) 

IMPLICIT  INTEGER-^  (I-N) 

REAL "8  AOS (3) ,AP AY (3) ,C (1100.0:5) ,A( 1100,0: 5 
INTEGER" 2  MOS f 1) , IZONeTI ) . ixfl) 

COMMON  CINF .BUDGET, TEST , NUMCEL 
READ  ( 10 , 100 )BUDGET 
READ  ( 10 , 200  jNUMCEL 
READ  ( 10 , 100)TIMAX 
READ  (  10 , 300) (AOS(I) ,1=1, 3) 

READ  (  10,300 1  ( APAY  ?I ) , I  =  1 , 3 ) 

READ  10, 100)BNSMAX 
READ  (10,400)PCTGFY 
READ  (10,400)0 
HLAM= - CINF 


5) , A ( 1 100 , 0 : 5 ) ,R(0:5) 


^  4  LJ  i.  •  J 

10 . 100 ) BUDGET 
10 , 200 1NUMCEL 

10 . 100) TIMAX 

10 . 300)  (AOS (I 

10 . 300)  (APAY(: 

10. 100 ) BNSMAX 

10 . 400 )  PCTGFY 

10 .400) Q 


1  =  1.3} 
,1  =  1,3 


SCALE= 1000 . DO 
BUDGET  =  BUDGET / SCALE 
BNSMAX=BNSMAX* PCTGFY 
TEST  =  CINF / 1 . IDO 
READ ( 1 1 , 500 ) ( MOS ( I ) , 


READ( 11,500) (MOS (I) , IZONE ( I ) , IX ( I ) , A (1,0), 1=1, NLD 
READ  ( 14,600 ) (A(I, 1) ,A(I ,2) ,1=1, NUMCEL ) 

DO  50  1=1, NUMCEL 


NUMCEL) 


MAXMUL=  6  - IZONE (I) 

READ  (15,700)(C(I,J) ,J=0, MAXMUL ) 
CONTINUE 

READ  ( 12 , 800 ) (A( I , 3 ) , 1= 1 .NUMCEL) 
READ  (13,900) (A(I, 4) ,1=1, NUMCEL) 

DO  150  1=1. NUMCEL 
MAXMUL=  6-1  ZONE ( I ) 

WT=A(I ,0) 

E=A(I,1) 

TE=A( 1,2) 

ACT= A ( 1,3) 


TI=A( I , 4 ) 

IF  (TI  .LT.  I. DO)  TI=50 . DO 
DO  7  5  J  =  0  , MAXMUL 
R(J)=C(I,J) 

CONTINUE 

AVGPAY= APAY ( IZONE ( I ) ) 

avgaos=aos?izoneTi) ) 

IF  (ACT  .GT.  0.0D0  THEN 
CW=WT*(TI/ UMAX)  /ACT 

ELSE 

CW=WT"TI/TIMAX 
END  IF 

IF  (E  .LT.  I. DO)  E= 1 . DO 
DO  125  J  =  0 , 5 

IF  (J  .LE.  MAXMUL)  THEN 
DlFTRM=TE-E"'R(  J ) 

IF  (DIFTRM  .LE.  O.DO)  THEN 

C  ( I ,  J  )  = CW"Q" (DIFTRM" "2 ) 

£  Xj  ^  ^ 

C ( I , J ) = CW* ( DIFTRM" "2 ) 

END  IF 

BONUS  =  PCTGFY " AVGPAY- AVGAOS " J 
IF  (BONUS  .GT.  BNSMAX )  BONUS^BNSMAX 
A(I .  J  )  =  rT  J  )'"E" BONUS  /  SCALE 
IF  (J  .GT.  0)  THEN 

HL=(C(I,0i-C(I.J))/A(I,J) 

IF  (HL  .GT.  HLAM)  HLAM=HL 
END  IF 
ELSE 

C(I , J)=CINF 
A  (I ,  J )  =  CINF 
END  IF 
CONTINUE 
CONTINUE 
FORMAT  (F12.0) 

FORMAT (14) 

FORMAT (3F7.0) 

FORMAT ( F4 . 2 ) 

FORMAT  (I5.2X.I1.2X, II, 2X,F4. 2) 

FORMAT (8X,2F10.0) 


RETURN 

END 

SUBROUTINE  BOUND ( ZL , IX , ZU , IXBEST , HLAM , C , A , COST , MOS ) 
IMPLICIT  REAL "8  (A-H.O-Z) 


'-/'-'i  HA'-'*.*  V_/ 

EPS = HLAM 7 1 . D05 
ENDL=0 .DO 


ENDL=0 .DO 
ENDR= . 0  IDO "HLAM 
ZUBEST  =  CINF 

AMBDA= ( ENDL  +  ENDR ) / 2 . DO 

CALL  MINFCN ( AMBDA . ZL , IX , ZU , COST , C . A , MOS ) 
IF (COST  .LE.  BUDGET)  THEN 
ENDR= AMBDA 

IF  (ZU  .LE.  ZUBEST)  THEN 
ZUBEST  =  ZU 
FLAM= AMBDA 
END  IF 
GOTO  100 
END  IF 
ENDL= AMBDA 
ENDR= 10.0  IDO "ENDR 
AMBDA= ( ENDL  +  ENDR ) / 2 . DO 

CALL  MINFCN ( AMBDA, ZL, IX , ZU , COST , C , A , MOS ) 


IF (COST  . LE.  BUDGET)  THEN 
ENDR= AMBDA 

IF  (ZU  .LE.  ZUBEST)  THEN 
ZUBEST = ZU 
FLAM= AMBDA 
END  IF 
GOTO  100 
END  IF 
ENDL= AMBDA 
ENDR=10 . 01D0--ENDR 


CALCULA 


AMBDA = ( ENDL+  ENDR ) / 2 . DO 

CALL  MINFCN( AMBDA , ZL , IX , ZU , COST , C , A , MOS ) 
IF (COST  .LE.  BUDGET)  THEN 
ENDR = AMBDA 

IF  (ZU  .LE.  ZUBEST)  THEN 
ZUBEST=ZU 
F LAM = AMBDA 
END  IF 
FT  SF 

ENDL= AMBDA 
END  IF 

IF  ((ENDR-ENDL)  . GT .  EPS)  GOTO 
ULATE  BEST  FEASIBLE  SOLUTION 


ULAtE  BEST  FEASIBLE  SOLUTION 

CALL  MINFEA ( FLAM , ZU , IXBEST , COST , C , A , IX , MOS ) 

RETURN 

END 

SUBROUTINE  MINFCN (AMBDA ,  ZL , IX , ZU , COST , C , A , MOS ) 
IMPLICIT  REAL--8  (A-H.O-Z) 

IMPLICIT  INTEGER-^  (I-N) 

REAL--8  C  f  1 100,0  :  lj,A(  1100,0:  1) 


REAL--8  C(1100,0:  1)  , A( 110 
INTEGER*^  IXliKMOS(I) 
COMMON  CINF, BUDGET, TEST, 
C0ST=0 . DO 
CELTOT  =  0 . DO 
ZU  =  0 . DO 

DO  400  1=1 .NUMCEL 
IF  (MOS(I)  .LT 
CELMIN=C ( 

COB J  =  C ( I  . 

INDEX= IX l 


NUMCEL 


N=C ( I , IX 


0)  THEN 

IX  j I ) ) + AMBDA" A ( I , IX ( I  )  ) 


GOTO  300 
END  IF 
CELMIN=CINF 
INDEX  =  0 
DO  200  J  =  0 , 5 

IF  (C(I.J)  .LE.  TEST) 
ZJ  =  C(I  ,  J  ) + AMBDA" A ( 
CELOBJ  =  C ( I , J | 

IF  (ZJ  .LE.  CELMIN 
IX(I)=J 
CELMIN=  ZJ 
COB J  =  CELOBJ 
INDEX= J 
END  IF 


THEN 

I,J) 


ELMIN )  THEN 


END  IF 
END  IF 
CONTINUE 

CELTOT  =  CELTOT  *  CELMIN 
ZU=ZU-COBJ 

COST=COST*A( I , INDEX) 

CONTINUE 

ZL= CELTOT -AMBDA "BUDGET 

RETURN 

END 

SUBROUTINE  MINFEA ( FLAM , ZU . IXBEST , COST , C , A , IX , MOS ) 
IMPLICIT  REAL" 8  (A-H.O-Z) 

IMPLICIT  INTEGER-4  (I-N) 

REAL" 8  C(  1100,0: 1 ) , A ? 1 100 , 0 : 1 ) 

INTEGER" 2  IXBEST ( 1) ,IX(  1) ,MOS( 1) 


COMMON  CINF , BUDGET .TEST , NUMCEL 
COSTCO. DO 
CELTOT  =  0 . DO 


THEN 

)) 


=  0  .  DO 

400  1=1. NUMCEL 
IF  (MOS(I)  .LT.  0)  THE 
COBJ=C (I . IX ( I )  ) 
INDEX= IX [ I ) 

ixbestTi J=IX(I) 
GOTO  300 
END  IF 
CELMIN=CINF 
INDEX =0 
DO  200  J  =  0 ,5 

IF  (C(I.J)  .LE.  TEST 
ZJ=Cf I , J)-FLAM*A( 
CELOBJ  =  C  ( I ,  J  ) 


IF  (ZJ  .LE.  CELMIN)  THEN 
IXBEST ( I )  =  J 
CELMIN=ZJ 
COB J  =  CELOBJ 
INDEX= J 
END  IF 
END  IF 
CONTINUE 
ZU  =  ZU  + COBJ 

COST  =  COST ♦ A ( I , INDEX) 

CONTINUE 

RETURN 

END 

SUBROUTINE  HEURIS (C .A, IXBEST , ZH . COST , MOS ) 
IMPLICIT  REAL" 8  (A-H.O-Z) 

IMPLICIT  INTEGER"4  (I-N) 

REAL  ••8  Cjf  1100.0:  1)  ,A(  UOO.O:  1)  ,UARRAY(  1100  ) 
INTEGER'-':?  IXBEST ft )  .MOS  (  1 ) 

COMMON  CINF, BUDGET, TEST, NUMCEL 
COSTl=COST 
RES ID = BUDGET- C0ST1 
DO  200  1  =  1  , NUMCEL 


■T .TEST .NUMCEL 


IF  ( MOS ( I )  .LT.  0) 

UARRaY (I)=-CINF 
GOTO  260 
END  IF 

IF  ( IXBEST ( I )  .GE.  5 
. GT.  TEST)  THEN 

UARRAY (I )= -CINF 
GOTO  2<J0 
END  IF 

RNUM=C(I . IXBEST (I ) ) -( 


THEN 

NF 


C ( I , ( IXBEST ( I ) 


RNUM=  C ( I . IXBEST (I ) ) - C ( I .  f I XBEST 
RDENOM  =  A[ I , ( IXBEST (I)* l) ) -  A ( I. I 
IF  (RDENOM  .GT.  RES  ID  .OR.  RNUM 
UARRAY ( I ) = -CINF 
GOTO  20O 

ELSE 

UARRAY ( I )=RNUM/ RDENOM 
END  IF 

IF  (UARR«\Y(I)  .LT.  0  .  DO  )  UARRAY 
CONTINUE 
UBEST= - CINF 
DO  300  1=1, NUMCEL 

IF  (UARRAY ( I )  .GT.  UBEST)  THEN 
UBE ST  =  UaRRAY ( I ) 

INDEX= I 
END  IF 
CONTINUE 

IF  (UBEST  .LT.  (-TEST)l  GOTO  400 
C0ST1  =  C0ST1-A( INDEX, I XBEST ( INDEX  )  ) 
"  +  A ( INDEX , ( IXBEST ( INDEX)  *  1  ) 


:  RNUM/ RDENOM 


I  XBEST (I ) ♦ 1 ) 
-All . IXBEST ( 

.  RNUM  . LE .  1 


. D -  0 8 )  THEN 


0 . DO )  UARRAY ( I ) = - UARRAY l I ) 


VRRAY  ( I ) 


UBEST)  THEN 


REENLISTMENT  RATE  FILE 


1833 

1 

0.10 

0.20 

0.21 

0.34 

.43 

1833 

2 

0 .65 

0.52 

0 . 69 

0.88 

.  9 

1833 

3 

0 .85 

0.90 

0.90 

0 . 90 

2111 

1 

0.24 

0.30 

0.36 

0.63 

0.67 

2111 

2 

0.50 

0 .51 

0.70 

0.80 

0.80 

2111 

3 

0.85 

0.90 

0.90 

0 . 90 

2112 

1 

0.24 

0 . 30 

0.36 

0 . 60 

0 .67 

2112 

2 

0.50 

0.51 

0 . 70 

0 . 80 

0.80 

2112 

3 

0 . 85 

0 . 90 

0.90 

0.90 

2131 

1 

0.24 

0.30 

0.63 

0 . 70 

0.80 

2131 

2 

0.50 

0.51 

0. 70 

0.80 

0.80 

2131 

3 

0.85 

0.90 

0.90 

0.90 

2142 

1 

0 . 20 

0.21 

0.24 

0 . 39 

0 .53 

2142 

2 

0.50 

0 .51 

0 . 70 

0 . 80 

0.80 

2142 

3 

0 .85 

0 . 90 

0.90 

0 . 90 

2144 

1 

0.20 

0.27 

0.36 

0 . 48 

0.56 

2144 

2 

0.50 

0.51 

0. 70 

0.80 

0.80 

2144 

3 

0.85 

0  90 

0.90 

0 . 90 

2145 

1 

0.20 

0.21 

0  .  37 

0 . 38 

0 . 65 

2145 

2 

0 . 50 

0.51 

0 . 70 

0.80 

0.80 

2145 

3 

0 .85 

0.90 

0.90 

0 . 90 

2146 

1 

0.24 

0.30 

0.36 

0.62 

0.65 

2146 

2 

0.50 

0.51 

0. 70 

0.80 

0.80 

2146 

3 

0 .85 

0 . 90 

0.90 

0 . 90 

2147 

1 

0 . 24 

0 . 30 

0  .  36 

.630 

.67 

2147 

2 

0 . 50 

0.51 

0 . 70 

0 .80 

0.80 

2147 

3 

0.85 

0.90 

0.90 

0.90 

2161 

1 

0.20 

0.30 

0.41 

0.42 

0.43 

2161 

2 

0 . 50 

0.51 

0. 70 

0.80 

0.80 

2161 

3 

0.85 

0 . 90 

0 . 90 

0 . 90 

2171 

1 

0.20 

0.30 

0.31 

0.32 

.43 

2171 

2 

0.50 

0.51 

0.70 

0.80 

0.80 

2171 

3 

0.85 

0.90 

0.90 

0.90 

2172 

1 

0.20 

0.27 

0.36 

0.48 

0.56 

2172 

2 

0.50 

0.51 

0 . 70 

0.80 

0.80 

2172 

3 

0.85 

0.90 

0.90 

0.90 

2311 

1 

0.23 

.25 

0.26 

0.41 

.  53 

2311 

2 

0.53 

0.70 

0.71 

0.72 

.8 

2311 

3 

0.85 

0 . 90 

0 . 90 

0.90 

2336 

1 

.21 

.22 

.  5 

0.94 

1 . 00 

2336 

2 

0.86 

.900 

.910 

.  920 

.  95 

2336 

3 

0.85 

0 . 90 

0.90 

0.90 

2512 

1 

0.21 

0.39 

0.40 

0 . 70 

.71 

2512 

2 

0.48 

.  75 

.  76 

.  88 

.  9 

2513 

1 

0.21 

0.39 

0.40 

0.70 

.71 

2513 

2 

.48 

.  75 

.  76 

.88 

.  9 

2519 

1 

0.21 

0 . 39 

0.40 

0 . 70 

.  71 

2519 

2 

0.77 

.  78 

.  8 

.  88 

.  9 

2519 

3 

0.85 

0 . 90 

0 . 90 

0.90 

2531 

1 

0.21 

0.39 

0.40 

0. 70 

.71 

.  5 
.  7 
.  7 
.85 
.  63 
0.63 
.  66 
.  66 
.  7 

0.63 

0.57 

0.63 

.  55 

1. 

.75 
.  75 
.75 


.75 


CELL  SIZE  FILE 


0121 

1 

1098 

0121 

2 

292 

0121 

3 

4 

0131 

1 

1151 

0131 

2 

301 

0131 

3 

11 

0151 

1 

3567 

0151 

2 

1036 

0151 

3 

102 

0161 

1 

194 

0161 

2 

104 

0161 

3 

38 

0193 

1 

14 

0193 

2 

422 

0193 

3 

640 

0211 

1 

18 

0211 

2 

57 

0211 

3 

24 

0231 

1 

256 

0231 

2 

221 

0231 

3 

59 

0241 

1 

24 

0241 

2 

66 

0241 

3 

30 

0251 

1 

24 

0251 

2 

76 

0251 

3 

21 

0311 

1 

13370 

0311 

2 

1010 

0311 

3 

34 

0313 

1 

252 

0313 

2 

62 

0313 

3 

2 

0331 

1 

2344 

0331 

2 

252 

0331 

3 

7 

0341 

1 

2600 

0341 

2 

238 

03 

41 

3 

6 

03 

51 

1 

2420 

03 

51 

2 

245 

03 

51 

3 

8 

03 

52 

1 

1034 

03 

52 

2 

252 

03 

52 

3 

5 

03 

69 

2 

1010 

03 

69 

3 

831 

04 

11 

1 

419 

04 

11 

2 

176 

TRAINING  COST  FILE 


0844 

2 

5211 

0844 

3 

5211 

0847 

1 

14027 

0847 

2 

14027 

0847 

3 

14027 

0848 

2 

15600 

0848 

3 

15600 

0861 

1 

16000 

0861 

2 

16000 

0861 

3 

16000 

1141 

1 

3256 

1141 

2 

3256 

1141 

3 

3256 

1142 

1 

9865 

1142 

2 

9865 

1142 

3 

9865 

1161 

1 

3484 

1161 

2 

3484 

1161 

3 

3484 

1169 

3 

12775 

1171 

1 

4386 

1171 

2 

4386 

1171 

3 

4386 

1181 

1 

3500 

1181 

2 

3500 

1181 

3 

3500 

1182 

1 

3500 

1182 

2 

3500 

1182 

3 

3500 

1316 

1 

7041 

1316 

2 

7041 

1316 

3 

7041 

1341 

1 

4536 

1341 

2 

4536 

1341 

3 

4536 

1345 

1 

8853 

1345 

2 

8853 

1345 

3 

8853 

1349 

3 

25216 

1371 

1 

2826 

1371 

2 

2826 

1371 

3 

2826 

1391 

1 

6705 

1391 

2 

6705 

APPENDIX  C 

SAMPLE  FROM  OUTPUT  FILE 

OBJECTIVE  FUNCTION  VALUE  IS:  1768.541 

THIS  SOLUTION  IS  WITHIN  0.00006802  70  OF  THE  OPTIMAL 

MOS  ZONE 

MULTIPLIER 

APPENDIX  D 
USER  INSTRUCTIONS 

A.  DATA  FILES 

Input  to  the  program  is  provided  through  six  data  files. 
File  names , specif ications  and  entries  are  detailed  in  the 
accompanying  tables.  In  general,  data  values  should  be 
right  justified  within  the  fields.  Rates,  percentages,  and 
the  exogenous  weighting  factors  may  be  placed  anywhere 
within  the  specified  field  provided  they  are  entered  as  two 
place  decimals  (e.g.,  .25,  .03  or  1.15).  With  the  exception 
of  the  Parameter  file,  all  files  contain  the  same  entries  in 
the  first  8  columns.  The  program  reads  data  values  sequen¬ 
tially  from  the  files.  Therefore,  it  is  vital  that  the 
MOSs  ,  Zones  and  accompanying  entries  from  one  file  corre¬ 
spond  line  for  line  with  all  other  files.  For  example,  if 
line  15  from  one  file  contains  entries  pertaining  to  MOS 
0193,  Zone  C,  then  line  15  from  all  other  files  should  refer 
to  MOS  0193,  Zone  C. 

Those  cells  for  which  a  non-zero  SRB  multiplier  is 
expressly  prohibited,  such  as  the  cells  in  the  Band  and 
Marine  Corps  Exchange  MOSs,  should  be  eliminated  from  the 
data  files.  Similarly,  there  are  MOSs  held  only  by  very 
senior  Marines  (E-8's  and  E-9's)  not  eligible  for  bonuses 
due  to  years - o f- service  requirements  exceeding  14  years. 
These,  too,  should  be  purged  from  the  data  files.  As 
written,  SRBMULT.FOR  can  read  a  maximum  of  1100  cells.  If  at 
some  point  it  becomes  necessary  to  increase  this,  it  will 
require  changing  the  array  and  matrix  dimensions  within  the 
program  from  1100  to  the  desired  number. 


.vnv.v.v  v''  w- 


TABLE  la 

DATA  FILE  SPECIFICATIONS 
Reenlistment  Plan  File  (RPLAN.DAT) 


Line 

Column 

Data  Description 

all 

1-5 

MOS 

6-7 

blank 

8 

Zone  (1  for  A,  2  for  B,  3  for 

9-18 

Number  of  Marines  eligible  for 
reenlistment  in  cell 

19-28 

Reenlistinent  target  for  cell 

TABLE  lb 

DATA  FILE  SPECIFICATIONS 
Cell  Size  File  (ACTNUM.DAT) 


Line  Column  Data  Description 


all 


6-7 

8 

9-14 


MOS 

blank 

Zone  (1  for  A,  2  for  B,  3  for  C) 
Total  number  of  Marines  in  cell 


TABLE  le 

DATA  FILE  SPECIFICATIONS 


Cell 

Weight  File  (WEIGHT . DAT ) 

Co lumn 

Data  Description 

1-5 

MOS 

6-7 

blank 

8 

Zone  (1  for  A.  2  for  B.  3  for  C) 

9-10 

blank 

11 

Preset  SRB  multiplier  (0  if 
preset) 

not 

12-13 

blank 

14-18 

Exogenous  weighting  factor  1 
ally,  use  l.Ou.  If  other  va; 

[  typic- 

Lues  a: 

used,  include  decimal,  e.g. 

.50) 

TABLE  If 

DATA  FILE  SPECIFICATIONS 
Parameter  File  (PARAM.DAT) 


Column  Data  Description 


1-12 


1-4 
1- 12 
1-7 
8-14 
15-21 


1-7 

8-14 

15-21 


1-12 

1-4 


1-4 


Budget  available  for  current  year 
bonuses  (dollars) 

Total  number  of  cells 
Maximum  training  cost  (dollars) 
Zone  A  avg  reenlistment  duration 
Zone  B  avg  reenlistment  duration 
Zone  C  avg  reenlistment  duration 
(years ) 

Zone  A  avg  monthly  base  pay 
Zone  B  avg  monthly  base  pay 
Zone  C  avg  monthly  base  pay 
(do  1 lars ) 

Maximum  permissible  bonus  (dollar 
Percentage  of  bonus  awarded  as  lui 
sum  (in  decimal  form,  e.g.,  .75) 

Over/under  factor 


B.  RUNNING  THE  PROGRAM 


To  use  the  program,  the  user  should  possess,  at  a 
minimum,  the  IBM  PC  configuration  detailed  in  Chapter  IV  and 
a  FORTRAN  77  compiler.  The  program  was  developed  and  tested 
using  the  Ryan-McFarland  RM/FORTRAN  compiler.  However,  the 
code  used  in  the  program  is  sufficiently  generic  that  other 
FORTRAN  77  compilers  should  be  acceptable. 

A  diskette  ("program  diskette")  containing  the  following 
files  will  be  provided  to  the  Enlisted  Plans  Section, 
Manpower  Department,  HQMC :  HOWTO . DOC ,  RUNCHK.BAT, 

FILCHK . FOR ,  FILCHK.EXE,  RUNSRB.BAT,  SRBMULT . FOR ,  and 
SRBMULT.EXE.  The  file  HOWTO.DOC  contains  a  copy  of  oper¬ 
ating  instructions.  Files  with  filetype  .BAT  are  executive 
routines,  those  with  filetype  .FOR  are  FORTRAN  source  code 
and  files  with  filetype  .EXE  are  executable  programs.  Use  of 
these  programs  is  detailed  in  the  following  sections.  Prior 
to  running  any  of  these  programs  certain  parameters  within 
the  DOS  configuration  file  CONFIG.SYS  must  be  changed.  If 
this  file  does  not  already  exist,  it  will  be  necessary  to 
create  it.  The  RM/FORTRAN  that  was  used  to  compile 
SRBMULT.EXE  and  FILCHK.EXE  requires  that  the  number  of  files 
that  can  be  opened  concurrently  be  increased  to  at  least 
ten.  Therefore,  the  following  line  must  be  added  to 
CONFIG . SYS : 

FILES  =  10 

Similarly,  the  number  of  disk  buffers  allocated  by  DOS  at 
system  startup  must  be  increased  to  ten.  Accordingly,  add: 

BUFFERS  =  10 


to  CONFIG.SYS.  Once  these  changes  have  been  made  to 
CONFIG.SYS,  reboot  DOS. 
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1 .  Checking  Data  Files  for  Conformity 

The  FORTRAN  program  FILCHK.EXE  is  provided  to  check 
the  files  for  MOS  and  Zone  conformity  prior  to  using 
SRBMULT .  FOR .  To  use  FILCHK.EXE  on  a  system  with  two  disk 
drives,  place  the  program  diskette  in  B  drive  and  the  disk¬ 
ette  containing  the  data  files  ("data  diskette")  in  A  drive. 
Enter  the  following  DOS  commands: 

Path  =  B : \ 


The  first  command  will  cause  DOS  to  check  the  B  drive  for 
executable  files.  The  second  sets  the  A  drive  as  the  default 
drive.  At  the  system  prompt,  enter  RUNCHK  to  execute  the 
program  RUNCHK.BAT.  This  executive  routine  will  set  neces¬ 
sary  file  definitions  and  then  execute  FILCHK.EXE.  On  a 
system  with  a  hard  disk,  place  all  files  in  the  same  direc¬ 
tory  on  the  C  drive  (hard  disk)  and  then  execute  RUNCHK.BAT 
as  before. 

The  FILCHK  program  uses  the  sequence  of  MOSs  and 
Zones  in  the  file  WEIGHT.DAT  as  a  reference  set  and  checks 
the  other  four  data  files  for  conformance  to  it.  If  a  non- 
conforming  entry  is  discovered,  the  applicable  file  defini¬ 
tion  number  and  line  number  are  written  to  the  screen.  The 
remaining  files  are  then  checked.  Once  detected  errors  are 
corrected,  RUNCHK.BAT  should  be  executed  again  until  all 
files  conform  with  the  reference  as  evidenced  by  no  error 
returns . 

2.  Running  SRBMULT 

To  run  SRBMULT.EXE  on  a  system  with  two  disk  drives, 
place  the  program  diskette  in  the  B  drive  and  the  data  disk¬ 
ette  in  the  A  drive.  Enter  the  following  DOS  commands: 

PATH  =  B:\ 

A: 

RUNSRB 


Several  minutes  later,  an  output  file  MULPLN.DAT  will  be 
created  on  the  diskette  in  the  A  drive.  It  will  contain  an 
MOS  and  Zone  listing,  the  recommended  SRB  multipliers,  the 
objective  function  value  and  the  percent  deviation  from  the 
optimal.  The  objective  function  value  is  useful  in  judging 
the  effect  of  introducing  exogenous  weighting  factors  or 
presetting  multipliers.  Instructions  on  how  to  do  both  of 
these  is  contained  in  the  following  section. 

Running  the  program  is  somewhat  simpler  on  a  system 
with  a  hard  disk.  All  files  should  be  placed  in  a  single 
directory  on  the  C  drive  (hard  disk)  and  that  drive  should 
be  made  the  default  drive.  SRBMULT.EXE  is  initiated  by 
typing  RUNSRB  to  execute  RUNSRB.BAT. 

C.  EXOGENOUS  WEIGHTING  FACTORS,  PRESET  MULTIPLIERS  AND  THE 

OVER/ UNDER  FACTOR 

Inclusion  of  an  exogenous  weighting  factor  (columns 
14-18  in  WEIGHT.DAT)  permits  the  user  to  exercise  additional 
control  in  weighting  the  effects  of  reenlistment  shortfalls 
in  particular  cells  based  on  his  judgement.  Weighting 
factors  W^  greater  than  the  default  value  of  1.0  result  in 
comparatively  greater  penalties  for  failing  to  meet  reen- 
listinent  targets.  Conversely,  weighting  factors  less  than 
1.0  result  in  lesser  penalties.  The  effect  of  the  factor  is 
linear.  Thus  a  2.0  results  in  a  penalty  twice  as  large  as  a 
1.0  and  a  factor  of  .5  reduces  the  penalty  by  one-half.  As 
a  weighting  scheme  based  on  training  costs  and  manning 
levels  is  already  utilized,  this  capabililty  should  be  used 
with  caution.  In  general,  a  value  of  1.0  should  be  used  for 
all  cells. 

The  user  can  also  preset  SRB  multipliers.  This  might  be 
necessary  if  factors  other  than  those  considered  by  the 
model  dictate  that  particular  cells  be  assigned  specific 
multipliers.  This  is  accomplished  by  entering  a  minus  sign 
(-)  in  column  1  of  the  pertinent  cell's  data  entry  in  the 


file  WEIGHT.DAT.  The  desired  multiplier  is  then  entered  in 
column  11  of  the  same  line.  In  the  output  file,  the  minus 
sign  will  appear  next  to  the  appropriate  MOS/Zone  thereby 
"flagging"  those  cells  with  predetermined  multipliers. 

The  last  entry  in  PARAM.DAT  is  the  over/under  factor  Q, 
the  weighting  factor  employed  to  permit  the  user  to  estab¬ 
lish  the  relative  importance  of  overshooting  or  under¬ 
shooting  reenlistment  targets.  For  example,  a  over/under 
factor  of  .70  implies  that  the  penalty  incurred  in  the  model 
for  overshooting  the  reenlistment  target  by  some  number  is 
only  70%  of  the  penalty  that  would  be  incurred  for  under¬ 
shooting  the  target  by  the  same  number.  Since  an  overage  in 
a  particular  cell  is  generally  considered  more  acceptable 
than  a  shortage,  the  factor  will  typically  be  less  than  1. 
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